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1. Introduction 

In this paper we are interested in a Levy process with no positive jumps, possibly indepen- 
dently killed at a constant rate, henceforth denoted by £ = 0) with law P. That is to 
say, under P, £ is a stochastic process which has almost surely cadlag paths, increments that are 
stationary and independent and killed at an independent rate k > 0. The case that k = cor- 
responds to no killing. Whilst it is normal to characterise Levy processes by their characteristic 
exponent, in the case that the jumps are non-positive one may also comfortably work with the 
Laplace exponent which satisfies, 

E(e**) = e^W, t > 0, 

where u > 0. Indeed, it is a well established fact that the latter Laplace exponent is strictly 
convex on [0, oo), see for example Bertoin pQ. Moreover, the Levy-Khintchine formula necessarily 
allows us to write tp in the form 

(1.1) if>(u) = -k + au + -a 2 u 2 + / (e M - 1 -«xl(i a; | < i))n(da;), 

for u > where k > 0, a € P, a 2 > and II is a measure concentrated on (— oo,0) satisfying 
/(-oo0)(l ^ x2 )n(d2;) < oo. Note in particular that our definition includes the case that — £ is 
a (possibly killed) subordinator. Indeed, when II satisfies j(l A |x|)II(dx) < oo and we 
choose o" = and a = — d + Q ^ xll(dx) we may write for u > 0, 



V>(u) = -K-du- (1 - e~ ux )u{dx), 

J(0,oo) 

where v(x, oo) = EE(— oo,— x). When d > 0, writing St = — £t for t > 0, d and v should be 
thought of as the drift and Levy measure of the subordinator S = (St,t > 0) respectively. 
Moreover, writing <j)(u) = —ip{u) for u > 0, we may think of (j) as the Laplace exponent of S in 
the classical sense, namely 

E(e~" St ) = e ~^ u)t , t > 0. 
l 



In general we shall refer to *$>i as the family of Laplace exponents of (possibly killed) Levy 
processes with no positive jumps which are killed at rate k > and are well defined on (£, oo) 
for I < 0. 

Our main objective is to introduce a parametric family of linear transformations which serves 
as a mapping from the space of Laplace exponents of Levy processes with no positive jumps 
into itself and therewith explore how a family of existing results for Levy processes may be 
extrapolated further. The paper is structured as follows. In the next section we introduce our 
three-parameters transformation and derive some basic properties. We also describe its connec- 
tion with some transformations which have already appeared in the literature. The remaining 
part of the paper deals with the applications of our transformation to different important issues 
arising in the framework of Levy processes and related processes. More specifically, in the third 
section we provide some ways for getting new expressions for the so-called scale function of Levy 
processes. Section 4 is devoted to the exponential functional of Levy processes and finally in the 
last section we develop some applications to the study of positive self-similar Markov processes. 

2. The transformation 7^/3 

We begin with the definition of our new transformation and consider its properties as a 
mapping on 

Definition 2.1. Suppose that ij) G where t < 0. Then for 5,(3 > 0, with the additional 
constraint that ifj'(Q+) = ip(0) = if f3 = 0, let 

Let us make some immediate observations on the above definition. Firstly note that Ts^ is a 
linear transform. In the special case that 5 = (3 we shall write 7/3 in place of 7/3,/? . The transform 
7/3 was considered recently for general spectrally negative Levy processes in Kyprianou and Patie 
|21j and for subordinators (as a result of a path transformation known as sliced splitting) in 
Gnedin [T7]. Next note that, for /3,7 such that (3 + 7 > 0, 

Note that one can in principle extend the definition of T$,p to Levy processes which have jumps 
on both sides which, for example, have all their exponential moments, by writing such processes 
as the difference of spectrally negative Levy processes and applying the 7^/3 linearly. 
In the special case that 5 = and /3 > we have £p := 7o,/3 satisfies 

£^(u) = i/>(u + 0)- V>(/3), u > -13, 

where, as usual, ip € This is the classical Esscher transform for Levy processes with no 
positive jumps expressed in terms of Laplace exponents. It will be convenient to note for later 
that if <l?(u) := ip(u)/u then we may write 

In particular we see that when j3 = 0, the assumption that ip'{ < d+) = allows us to talk safely 
about $(0+). 

One may think of 7^/3 as one of the many possible generalisations of the Esscher transform. 
For f3 > 0, the latter is a well-known linear transformation which maps into itself and has 
proved to be a very effective tool in analysing many different fluctuation identities for Levy 
processes with no positive jumps. It is natural to ask if 7j i( g is equally useful in this respect. A 
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first step in answering this question is to first prove that Tgp also maps ^ into itself. This has 
already been done for the specific family of transformations Tp in Lemma 2.1 of |21j . 



Proposition 2.2. Suppose that if) G ^fi where £ < 0. Fix S, ft > u>ii/i i/ie additional constraint 
that ip'(0+) = ip(0) = if (3 = 0. Then Tg pip € C ^ and has no killing component. 

Moreover, if tp has Gaussian coefficient a and jump measure U then Tg t pip also has Gaussian 
coefficient a and its Levy measure is given by 

e^ x U(dx) + 6e^ x U(x)dx + 5^e? x dx on (-oo, 0), 

P 

where H(x) = IT(— oo, —x) and we understand the final term to be zero whenever k = 0. 

Proof. Recall from earlier that Tg,pip(u) = £pip{u) — 5£p(p{u). Moreover, from its definition, it 
is clear that argument of Tg p may be taken for all u such that u + B > i. It is well understood 
that £pip is the Laplace exponent of a spectrally negative Levy process without killing whose 
Gaussian coefficient remains unchanged but whose Levy measure is transformed from II (dx) to 
e P x Yl(dx). See for example Chapter 8 of |20j. The proof thus boils down to understanding the 
contribution from —5£p(p(u). 

A lengthy but straightforward computation based on integration by parts shows that 

<K(n) = -- + (a-n(-l)) + Vn+ f (1 (N<1) - e xu )TL{x)dx. 

From this it follows that 

) = -S^-^—-6(a-U(l))u-5 f {1 - e xu )e^ x TL{x)dx 

/0 K roo 

(1 - e xu )e f3x Il{x)dx -5- (1 - e~ ux )e- px dx. 
-oo P JO 

Here we understand the final integral above to be zero if k = 0. In that case we see that 
—5£pip{u) is the Laplace exponent of a spectrally negative Levy process which has no Gaussian 
component and a jump component which is that of a negative subordinator with jump measure 
given by 5e^ x H(x)dx + (3~ 1 K5e l ^ x dx on (— oo, 0). □ 

Remark 2.3. As Tg p maps ^£ in to itself, it is possible to iterate the application of this transfor- 
mation to produce yet further examples of spectrally negative Levy process exponents. Consider 
for example the transformation Tgp = o Ts,p where 5, (3 > 0, 7 > 0. It is straightforward to 
check for this particular example that 

' p u + 7 \ u + y + p p J 

Whilst it is now clear that the mappings may serve as a way of generating new examples 
of Levy processes with no positive jumps from existing ones, our interest is largely motivated 
by how the aforesaid transformation interacts with certain path transformations and fluctuation 
identities associated to Levy processes. Indeed, as alluded to above, starting with Urbanik [35] . 
the formalisation of these transformations is motivated by the appearance of particular examples 
in a number of such contexts. On account of the diversity of these examples, it is worth recalling 
them here briefly for interest. 



Applications to pssMp: In Kyprianou and Patie [2JJ the transformation Tp turns out to give 
a natural encoding for the Ciesielski- Taylor identity as seen from the general point of view of 
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spectrally negative positive self similar Markov processes (pssMp). Recall that Lamperti |24| 
showed that for any x > 0, there exists a one to one mapping between Pi og: r> the law of a Levy 
process killed at rate k > 0, starting from logx, and the law Q x of a a-self-similar positive 
Markov process X on (0, oo) where a€M, i.e. a Feller process which satisfies for any c > 0, 

(2.1) ({X t ,t > 0), Q cx ) ( = } ((cX r . t) t > 0), Q x .) . 

More specifically, Lamperti proved that X can be constructed from £ as follows 

(2.2) log(X t ) = U> *>0, 
where 

A t = inf{s > 0; S s := / e a ^ du > t}. 
Jo 

In the case that £ is spectrally negative with Laplace exponent ip, we shall denote the law of 
X by Qx • Under reasonably broad conditions it is well understood how to extend the definition 
of X as Markov process to include into its state space; that is to say, how to give meaning to 
the law Qo- This is done either by constructing an entrance law for X at 0, or by constructing 
a recurrent extension from 0; see for example [H [I0l [32l I15| . Within the class of pssMps which 
are spectrally negative, the construction of a unique entrance law is always possible whenever 
V/(0+) > and k = 0. The existence of a unique recurrent extension when ip'(0+) < or k > 
requires that 6 := sup{A > : ^(A) = 0} € (0, a). 

Ciesielski- Taylor identity: Within this class of spectrally negative pssMps and in the spirit of 
the original Ciesielski- Taylor identity for Bessel processes (which themselves are included in the 
aforementioned category), [21] showed that for each a > 0, 

(T a ,Qt)= (jT ' l { x 3 <a}ds,®fA , 

where T a = inf{i > : Xf = a}. 

Continuous state branching processes with immigration: Here, we take X to be a self-similar 
continuous state branching process with immigration. We refer to Patie \27\ Section 4] for more 
detailed information on this family of Markov processes. In particular, it is shown that X is 
self-similar of index < a < 1 and its branching mechanism is tp{u) = — ^u a+1 , c > 0, while its 
immigration mechanism is given by x (u) = c5{^-)u a , 5 > 0. The Levy process £ associated to 
X via the Lamperti mapping (|2.2p is characterized by the Laplace exponent if)( a '^ of — £ which 
admits the following form 

where tp(u) = c j?fe^ is the Laplace exponent of the spectrally negative Levy process associated, 
via the Lamperti mapping, to the spectrally negative (a + 1)— stable Levy process killed upon 
entering into the negative real line, see [271 Proposition 4.11]. 

Sliced splitting of subordinators: In combination with the study of certain combinatorial 
structures, Gnedin [17] introduced the following method of sliced splitting the path of subordi- 
nators which corresponds to the application of a special case of the transformation we introduce 
here. 

For a subordinator (no killing) with Laplace exponent (f>, run its path until it first crosses a 
level e, which is independent and exponentially distributed with parameter (3. Remove the path 

4 



of the subordinator which lies above this level and issue an independent copy of the original sub- 
ordinator from the crossing time at the spatial position e. Repeat this procedure independently 
within the space-time framework of the latter subordinator and so on. 

To describe more precisely the resulting process, suppose that {ej : i > 1} are the inter- 
arrival times of arrival of a Poisson process with rate j3 and suppose that {S^ : i > 1} arc 
independent copies of the subordinator whose Laplace exponent is <j>. Let tq = and for i > 1, 
Tj = M{t > : sf ) > ej}. The path which was earlier described in a heuristic way is defined 
iteratively by So = 0, and, for i > 1, 

St = S Ti _ 1 + m.m{S^} Ti _ i ,e i } when t € (Ti_i,Tj], 

where Tq = and, for i > 1, Tj = t\ + • • • + t,. Gnedin [T7j shows that the resulting process 
5 = (St,t > 0) is again a subordinator and moreover its Laplace exponent, <f>p, satisfies 

for u > 0. 



3. Scale functions for spectrally negative Levy processes 

Scale functions have occupied a central role in the theory of spectrally negative Levy processes 
over the last ten years. They appear naturally in virtually all fluctuation identities of the latter 
class and consequently have also been instrumental in solving a number of problems from within 
classical applied probability. See Kyprianou [20J for an account of some of these applications. 
Despite the fundamental nature of scale functions in these settings, until recently very few 
explicit examples of scale functions have been found. However in the recent work of Hubalek 
and Kyprianou |18j . Chaumont et al. |14| . Patie |27| . Kyprianou and Rivero [23], many new 
examples as well as general methods for constructing explicit examples have been uncovered. 
We add to this list of contemporary literature by showing that the some of transformations 
introduced in this paper can be used to construct new families of scale functions from existing 
examples. 

Henceforth we shall assume that the underlying Levy process, £, is spectrally negative, but 
does not have monotone paths. Moreover, we allow, as above, the case of independent killing 
at rate k > 0. For a given spectrally negative Levy process with Laplace exponent ip, its scale 
function : [0, oo) h-» [0, oo) is the unique continuous positive increasing function characterized 
by its Laplace transform as follows. For any n > and u > 8 := sup{A > : ip(X) = 0}, 

[°° e- ux W^{x)dx = ^—. 
Jo VW 

In the case that k > 0, is also known as the K-scale function. 

Below we show how our new transformation generates new examples of scale functions from 
old ones; first in the form of a theorem and then with some examples. 

Theorem 3.1. Let x,/3>0. Then, 

(3.1) W T M = e- px W^{x) + (3 I e-^W^{y)dy. 

Jo 

Moreover, if ip'(0+) < 0, then for any x,/3,5 > 0, we have 

WrsM x ) = e- ex (w^(x) + 8e Sx J\- 5 yW^y)dy^ 
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Proof. The first assertion is proved by observing that 



OO 

ux 



e- ux W T ^(x)dx 



o ' u^{u + P) 

, ^ 1 P 

( 3 - 2 ) = T7 m + 



i/>[u+p) uip{u + py 

which agrees with the Laplace transform of the right hand side of (13. ip for which an integration 
by parts is necessary. As scale functions are right continuous, the result follows by the uniqueness 
of Laplace transforms. 

For the second claim, first note that Ts^tp = (u + 9 — 5)tp(u + 6)/(u + 9). A straightforward 
calculation shows that for all u + 5 > 9, we have 

( w u + 6 

w (x)cte = 

' uip{u + d) 

The result now follows from the first part of the theorem. □ 



(3.3) / e- ux e {e -^ x W Ts0 Jx)dx 



When ip(0+) > and ip(0) = 0, the first identity in the above theorem contains part of 
the conclusion in Lemma 2 of Kyprianou and Rivero [23J. However, unlike the aforementioned 
result, there are no further restrictions on the underlying Levy processes and the expression on 
the right hand side is written directly in terms of the scale function as opposed to elements 
related to the Levy triple of the underlying descending ladder height process of £. 

Note also that in the case that ip is the Laplace exponent of an unbounded variation spectrally 
negative Levy process, it is known that scale functions are almost everywhere differentiable and 
moreover that they are equal to zero at zero; cf. Chapter 8 of [20]. One may thus integrate 
by parts the expressions in the theorem above and obtain the following slightly more compact 
forms, 

W r ^(x) = / e-^W^(y)dy and W w (x) = e^" 5 )* / e^W^dy. 
Jo Jo 

We conclude this section by giving some examples. 

Example 3.2 ((Tempered) Stable Processes). Let ip KC (u) = (u + c) a — c a — k where 1 < a < 2 
and k, c > 0. This is the Laplace exponent of an unbounded variation tempered stable spectrally 
negative Levy process £ killed at an independent and exponentially distributed time with rate k. 
In the case that c = 0, the underlying Levy process is just a regular spectrally negative a-stable 
Levy process. In that case it is known that 

poo 



f°° 1 1 

JO U a — K 



and hence the scale function is given by 

W* Kfl (x) = i a " 1 ^, a (ra a ), x>0, 

where £ a ^(x) = Y^=o v(an+p) s t &n ds for the generalized Mittag-Leffler function. (Note in 
particular that when k = the expression for the scale function simplifies to T(a) x a ~ l ). 
Since 

/•OO 1 

it follows that 

WW» = e— W^ K+ca>0 {x) = e-^x^S^dK + c a )x a ) 



Appealing to the first part of Theorem 13. II we now know that for /? > 0, 

Jo 

Note that t/)' K c (0+) = ac Q_1 which is zero if and only if c = 0. We may use the second and 
third part of Theorem 13, II in this case. Hence, for any 5 > 0, the scale function of the spectrally 
negative Levy process with Laplace exponent Ts,otpo,o is 

wv s aoW = T^zrr/* £ e ~ Sy y a ~ 2d y 

<V ' -e 5x T(a- l,Sx) 



r«-i 



where we have used the recurrence relation for the Gamma function and T(a, b) stands for the 
incomplete Gamma function of parameters a, b > 0. Moreover, we have, for any /3 > 0, 

Finally, the scale function of the spectrally negative Levy process with Laplace exponent Tf ip K ,o 
is given by 



W Ts, ^, (x) = (x/5) a 1 £ a , a -i (x; j) 



where we have used the notation 

T(x;an + (3)K n 



£ a ,/3 {x; k) = 2^ — 



n=0 + 

Example 3.3 (T-stable process). Let us consider, for 1 < a < 2, tp(u) = (u — l) a , u > 0, where 
(u) a = -feffi stands for the Pochhammer symbol, and note that = 0. It was shown in 

[HI [27] that its scale function is given by 

W i ,(x) = -^-e- x (l-e- x ) a -\ x>0. 
Via) 



Hence, the scale function of the spectrally negative Levy process with Laplace exponent 7s, iV' is 

2 05-l)a 



WtsM x ) = Vr- / e-^y( e y - i) a - 2 ( a - e y)d y 







and 



w t s ^= = a + 1 _ g /* e-^ +1+ ^y( e y - i) a -\a - e y)d y 







p + l-o J 



4. Exponential functional and length-biased distribution 



In this part, we aim to study the effect of the transformation to the law of the exponen- 
tial functional of some Levy processes, namely for subordinators and spectrally negative Levy 
processes. We recall that this random variable is defined by 

/>oo 

I^p = e~^ s ds. 
Jo 

Note that lim^oo £ t = +00 a.s. 44> 1^ < 00 a.s. which is equivalent, in the spectrally negative 
case, to E[£i] = ifi'(0 + ) > 0. We refer to the survey paper of Bertoin and Yor [7] for further 
discussion on this random variable. We also mention that Patie in [28], [31] and [30J, provides 
some explicit characterizations of its law in the case £ is a spectrally positive Levy process. In 
the sequel, we denote by (dx) the law of 1^ . 



4.1. The case of subordinators. Let us first assume that £ is a subordinator, that is a non- 
negative valued Levy process. We recall that, in this case, the Laplace exponent of £ is given 
by 4>(u) = —tp(u), u > 0, where in particular we allow for the case of killing (k > 0). In that 
case we write I a, in place of 1^ and if 0(0) = k, we have that 

e K 

e~^ s ds 

where e K stands, throughout, for an exponential random variable of mean > 0, independent 
of £ (we have eo = 00). Before stating our result, we recall that Carmona et al. [13] determine 
the law of Ia through its positive entire moments as follows 

(4.1) m = m - Wyn = , 1 .... 

Theorem 4.1. For any k,/3>0, the following identity 

x Pf^j dx 
E[J; 

holds. 



(4.2) f T ^(dx) = *>0, 



Proof. Carmona et al. [13], see also Maulik and Zwart [261 Lemma 2.1], determine the law 
of 1^ by computing its positive entire moments which they derive from the following recursive 
equation, for any s, ft > and k > 0, 

(4-3) E[If^] = I^1e[I^} 



s + 13 

On the other hand, we also have, for any s,/3 > 0, 



We get the first assertion by uniqueness of the solution of the recursive equation and the fact 
that the law of 1a, is moment determinate. □ 



Before providing some new examples, we note from Theorem 14.11 that if Is = AB for some 
independent random variables A, B then the positive entire moments of Ifpfa /3 > admit the 
following expression 

(44) MI" 1 - E[An+ " ] E[Bn+ " ] n-01 

Example 4.2 (Poisson process). Let £ be a Poisson process with mean c = — log(g) > with 
< q < 1, i.e. 4>{u) = — log(g)(l — e~ u ),u > 0. Biane et al. [3] computed the law of 1$ by means 
of q— calculus. More precisely, they show that its law is self-decomposable and is given by 

OO n(n-l) 

U(dy) = £(-l)"e-^" \ ' dy, y > 0, 

where 

n— 1 oo 

(a;q)n = Y[{1 ~ aq J ), {q;q)oc = ~ a ^ J ) 

k=l k=l 

and its Mellin transform is given, for any s > 0, by 



The image of £ by the mapping 7g is simply a compound Poisson process with parameter 
c and jumps which are exponentially distributed on (0, 1) with parameter /?, i.e. Tp(j)(u) = 
log(g)(l — e~( u+ ^)). Thus, we obtain that the law of Ij-^ has an absolute continuous 
density given by 

£j OO n(n-l) 

f Tg jdx) = -g_ yv*/g"(-i)" q ^ i — -dx, x >o. 

To conclude this example, we mention that Bertoin et al. [3] show that the distribution of the 
random variable defined, for any bounded Borel function /, by 

with an independent copy of 1^, shares the same moments than the log normal distribution. 
It is not difficult to check that such transformation applied to It b 4> does not yield the same 
properties. 

Example 4.3 (Killed compound Poisson process with exponential jumps). Let £ be a compound 
Poisson process of parameter c > with exponential jumps of mean b^ 1 > and killed at a rate 
k > 0. Its Laplace exponent has the form (j)(u) = c-^^ + k and its Levy measure is given by 
v{dr) = cbe~ br dr, r > 0. We obtain from (|4.ip that 

E[jni n!r(n + fe+l)r( Kfe + l) 



{{K + c)) n T{b+l)T{n + K b + l) 
where we have set = -^h,b. Then, noting that b — Kb > 0, we get the identity in distribution 

I<i> - ((« + c))~ x exB{n h + 1,6- k 6 ) 

where -B(a, b) stands for a Beta random variable of parameters a, b > and the random variables 
on the right hand side are considered to be independent. The case k = was considered by 
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Carmona et al. [12]. Finally, observing that, for any (3 > 0, Tg0(u) = c u+ ^ +j3 + and its 

Levy measure is vp(dr) = e~^ r (c(b + (3)e~ br + n/3)dr, r > 0, we deduce from Theorem 14, 1| the 
identity 

I T? <p = (a(« + l)) _1 G(/3 + 1)B(k 6 + /3 + 1, 6 - « 6 ) 
where G(a) is an independent Gamma random variable with parameter a > 0. 

Example 4.4 (The a-stable subordinator) . Let us consider, for < a < 1, (j)(u) = u a , u > 
and in this case u(dr) = ^jjz^'dr,r > 0. The law of 1^ has been characterized by Carmona 
et al. [12, Example E and Proposition 3.4]. More precisely, they show that the random variable 
Z = log is self-decomposable and admits the following Levy-Khintchine representation 

logE[e iuZ ] = log(r(l + iu)) 1 - a 

r0 „s 



= (1 — a) I -iuEy + / (e ms - ius - 1)— —ds 

\ J—oo I s I (.1 — e j 

where E y denotes the Euler constant. First, note that Tp4>(u) = (u + /3) Q — /3(u + /S)"^ 1 and 
vp(dr) = Y(i-a) e~^ r r~(" +1 ^ (a + (3r)dr,r > 0. Thus, writing Z' 3 = log(/-7- 3 0), we obtain 

logE[e**'] = log^ (1 + /3 + mr "" 



r(i + /?) 

= (l-a)liuT(l+0) + J (e ms -ius-1) ds\ 

where T stands for the digamma function, T(z) = j&Q- Observing that lim a _>Q Tp(j)(u) = 
by passing to the limit in the previous identity we recover the previous example. 

Example 4.5 (The Lamperti-stable subordinator). Now let 4>(u) = (fto( u ) = (otu) a ,u > 0, with 
< a < 1. This example is treated by Bertoin and Yor in [TJ. They obtain 

r W 

1$ = e i e 1 . 
Hence, with Tp(f>(u) = ^(a(u + (3)) a , we get 

4.2. The spectrally negative case. Let us assume now that £ is a spectrally negative Levy 
process. We recall that if E[£i] < 0, then there exits 9 > such that 

(4.5) E[e^] = 1 

and we write ipg(u) = ip(u + 8). We proceed by mentioning that, in this setting, Bertoin and Yor 
[6] determined the law of /</, by computing its negative entire moments as follows. If ip(0) = 
and ip'(0 + ) > 0, then, for any integer n > 1, 



(4.6) E[/-"] = V'(0 + ) 

with the convention that Ef/^" 1 ] = t//(0 + 
for Levy processes yields, for any a > 



with the convention that Ef/^" 1 ] = t//(0 + ). Next, it is easily seen that the strong Markov property 



1^ ( J e -t°ds + e- a lL 
Jo 
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where T a = inf{s > 0; £ s > a} and IL is an independent copy of 1^, see e.g. [32]. Consequently, 
1$ is a positive self-decomposable random variable and thus its law is absolutely continuous 
with an unimodal density, see Sato [33]. We still denote its density by f^p(x). Before stating our 
next result, we introduce the so-called Erdelyi-Kober operator of the first kind and we refer to 
the monograph of Kilbas et al. [19] for background on fractional operators. It is defined, for a 
smooth function /, by 

B a > 5 f(x) = [ X r a (x-r) S - l f(r)dr, x > 0, 

T(o) Jo 

where £He(<5) > and 9\z{a) > 0. Note that this operator can be expressed in terms of the 
Markov kernel associated to a Beta random variable. Indeed, after performing a change of 
variable, we obtain 

1 [a + o + 1) 

which motivates the following notation 

(4.7) D a ' 5 f(x) = r( " + ^t, 1) B^/(x). 

T[a + 1) 

Theorem 4.6. (1) If tp'(0 + ) > 0, then for any (3 > 0, we have 

/ r ,*(s) = ^^,x>0. 

In particular, 1^ is the length-biased distribution of I-r^. 
(2) Assume that ip'{Q + ) < 0. Then, for any < 5 < 0, we have 

(4.9) I Ts<Bi> ^ B-\9- 5,5)1^. 

This identity reads in terms of the Erdelyi-Kober operator as follows 

f Tsjj1 p(x) = V e - 5 - l > 5 U e (x), x>0. 
In particular, we have 



(4.10) f Ts ^{x) ~ r( ^_ § f K; S ] * M os x^oo, 



(f(t) ~ g(t) as t — > a means that lim^a = 1 for any a £ [0,oo]). Combining the 
two previous results, we obtain, for any < 5 < 6 and (3 > 0, 

X ~P 

fr? J x ) = , -a 
'jjr w It p 



and 



,. . . T{B) s-e-p-i 

. (x) ~ — - — ', n Fir x p as X ->■ CO. 



-Pi 



Remark 4.7. We point out, that under the Cramer condition, Rivero |32t Lemma 4] obtains 
the following large asymptotic for the tail of the distribution 

(4.11) P (h > x) ~ Cx~ e as x — y co, 

where C > and ^(6) = with > and ^ being here the Laplace exponent of a non-lattice 
Levy process. Rivero has a characterization of the constant C only in the case 9 < 1. In our 

li 



setting, it is plain that, for any j3 > 0, Tf e ip(S — 9 — (5) =0. Moreover, f^p ^(x) is unimodal 

and hence ultimately monotone at infinity, therefore our asymptotic results, up to the constant 
characterization, could have been also deduced from (|4.1ip . 

Proof. We start by recalling the following identity due to Bertoin and Yor [2] 

T it ( d ) -M 
mi<j> = e ei 

where <j)(u) = ip(u)/u, u > 0, that is cj) is the Laplace exponent of the ladder height process of the 
dual Levy process, M = sup t>0 {— ^} is the overall maximum of the dual Levy process and the 
random variables are considered to be independent. Thus, recalling that E[e _sM ] = ip'(0 + )/(j)(s), 
we have, for any s > 0, 

(s+ i) y/(o+)r( g + 2) 

^(0+)r(a + l) 

(4.12) = ^E[/-] 

where we have used the recurrence relationship satisfied by 1^, see (|4.3|) . Similarly to the case 
of subordinators, we have, for any s > 0, 



- *< 5 + «-E[/r 



and 

L 4, J g + p L v- 

The first claim follows. Next, we have both ^(0 + ) > and Ts,e^'(0 + ) = ^^(0 + ) > as 
5 < 6. Thus, the random variables and lT ss tp are weu defined. Moreover, from (|4.6p . we get, 
for any integer n > 1, 

-s) mz{Ts,em 

i^W - *) IDS 



T(n) 

r(n + 0- 5)T{6) UkZlMk) 



Y(9-5)Y(n + 8) T(n) 
The identity (|4.9p follows by moments identification. Then, we use this identity to get, for any 



x > 0, 



Y(6)Y(6 - 5) 

x- e T{9) 
Y{8)Y{6 - 5) 



,6»— 5-1 /_ „.\5-l 



(x-u) f^ e (u)du 
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Next, we deduce readily from (|4,12p that the mapping s *— > IE[^ e s ] is analytic in the right-half 
plane 9\e(s) > —6. In particular, for any < 5 < 9, we have Ef/^ 5 ] < oo. Then, the large 



asymptotic behavior of the density is obtained by observing that 

x' 



HM X ) = rmwl r( ^ / u»->{l-u/x)^f^(u)du 







- & - e - l T{6) 



u e s f 1 p e (u)du as x — >• oo, 



T(5)T{9 - 5) 

which completes the proof. □ 

Example 4.8 (the spectrally negative Lamperti-stable process). Let us consider the Lamperti- 
stable process, i.e., for 1 < a < 2, ip(u) = ((a — l)(u — 1)) Q , u > 0. Recall that = 0, 
tpi(u) = ((a — l)u) a and that this example is investigated by Patie [27]. We get 



E[/T n ] = no" 



(o-l) 



' ! ' ' ' V(a + 1 
Thus, 1^ = e 1 1 " ±; . Then, for any < 5 < 1, 



and for any (3 > 

( i fi(l+/3-,5, ( 5)- 1 G(/3 + l)-( Q - 1 ). 



5. Entrance laws and intertwining relations of pssMp 

In this part, we show that the transformations Ts t p appear in the study of the entrance law 
of pssMps. Moreover, also they prove to be useful for the elaboration of intertwining relations 
between the semigroups of spectrally negative pssMps. We recall that the Markov kernel A 
associated to the positive random variable V is the multiplicative kernel defined, for a bounded 
Borelian function /, by 

Af(x)=E[f(Vx)}. 

Then, we say that two Markov semigroups P± and Qt are intertwined by the Markov kernel A if 
(5.1) P t A = AQ t , t > 0. 

We refer to the two papers of Carmona et al. [TT] and |13] for a very nice account of inter- 
twining relationships. In particular, they show, by means of the Beta-Gamma algebra, that the 
semigroup of Bessel processes and the one of the so-called self-similar saw tooth processes are 
intertwined by the Gamma kernel. Below, we provide alternative examples of such relations for 
a large class of pssMps with stability index 1. Recall that the latter processes were defined in 
Section [2j We also mention that for any a > —1,5 > 0, the linear operator D a,i defined in 
(14. 7p is an instance of a Markov kernel which is, in this case, associated to the Beta random 
variable B(a + 1, 5). In what follows, when X is associated through the Lamperti mapping to a 
spectrally negative Levy process with Laplace exponent ip, we denote by pf its corresponding 
semigroup. When min(^(0), ip'(0 + )) < and 6 < 1, then pf stands for the semigroup of the 
unique recurrent extension leaving the boundary point continuously. Using the self-similarity 
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property of X, we introduce the positive random variable denned, for any bounded borelian 
function /, by 

P^/(0)=E[/(t^)]. 

Recall that Bertoin and Yor [5] showed that, when ip'(0 + ) > 0, the random variable J w is 
moment-determinate with 



(5.2) 



T(n+1) 



n 



1,2... 



Before stating the new intertwining relations, we provide some further information concerning 
the entrance law of pssMps. In particular, we show that the expression (|5.2p of the integer 
moments still holds for the entrance law of the unique continuous recurrent extension, i.e. when 
min (ip(0), ip'(0 + )) < with 6 < 1. We emphasize that we consider both cases when the process 
X reaches continuously and by a jump. 

Proposition 5.1. Let us assume that min (V'(O), il>'(0 + )) < with < 1. Then, we have the 
following identity in distribution 

(5.3) V= } B(l-e,0)/I Tl _^ e 

where B(a,b) is taken independent of the random variable Lj- 1 _ g ^ g . Moreover, the entrance law 
of the unique recurrent extension which leaves continuously a.s. is determined by its positive 
entire moments as follows 



T(n+1) ' 



n 



1,2... 



Proof. We start by recalling that Rivero [32 j, Proposition 3] showed that the (/-potential of the 
entrance law of the continuous recurrent extension is given, for a continuous function /, by 



,-6» /-oo 



e-9*E [f(tJ^)\ dt= H — 

L>8 Jo 



f(u)E 



u e du 



where we have used the self-similarity property of X and set Cg = T(l — 9)E 



. Performing 



the change of variable t = uL^ on the right hand side of the previous identity, one gets 



,-6» /-oo 



e-«*E [f(tJ^] dt= H — 

<^e Jo 



- qt E 



We ' We 



t~ y dt. 



Choosing f(x) = x s for some s £ iM. the imaginary line, and using again the self-similarity 
property, we get 

fOO 

'-'-to 



poo 

/ e -^E[j;] = r( s + i)^- 1 E[j^] 

Jo 



where we have used the integral representation of the Gamma function r(z) = J °° e t t z dt, < 0\t{z) > 
— 1. Moreover, by performing a change of variable, we obtain 

;>oo 

(5.4) J t s E[e- qtI ^]t- 1 - e dt = q- s - - 1 T(s-e + l)E[I- g s - e - 1 }. 

Putting the pieces together, we deduce that 

r(«-0 + i) 



(5.5) 



E[4] 



E\I, 



r(* + i)r(i-f) e[i 
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and the proof of the first claim is completed by moments identification. Next, we have 



E[J£] 



r(n + i 



■E 



/•- ra+0-1 



r(n + i)r(i - e)E 

r(n + i-g) 
r(n + i)r(i-e) 

T(n + 1 - 0) V(l) Ilfei 1 srf^C* + 1) 



r-n 



r(n + i)r(i — 6») 1 — 6» 
r(n + 1 - 0)r(n)r(2 



r(n) 

^(i) n£M* + i) 



r(n + l-0)r(n 



lira 



i 



r n 



61) and 



r(n + i) 

where we have used, from the second identity, successively the identities (|5.5[) 
the recurrence relation of the gamma function. We point out that under the condition 9 < 1, 
ip(k) > for any integer k > 1. The proof of the Proposition is then completed. □ 

Before stating our next result, we recall a criteria given by Carmona et al. [11^ Proposition 
3.2] for establishing intertwining relations between pssMps. If / and g are functions of Co(R + ), 
the space of continuous functions vanishing at infinity, satisfying the condition: 

(5.6) Vt>0, P t /(0) = P t g(0) then f = g. 

Then the identity (|5.ip is equivalent to the assertion, for all / € Co(R + ), 

PiA/(0) = Qi/(0). 
Finally, we introduce the following notation, for any s € C, 

M^s) = E[4 s ). 

Theorem 5.2. (1) Assume that -0'(O + ) < and M^ e (s) ^ for any s € R. Then, for any 
5 < 9 + 1, we have the following intertwining relationship 

t > 

and the following factorization 

(d) 



(5.7) 



B(l + 0-6, 5)J^ e 



holds. 



(2) Assume that min (ip(0), i>'(0 + )) < with 9 < 1 and Mr_ e ^ g {s) / for any s G R. Then, 
we have the following intertwining relationship 



1,0 T}4> 



D ■ P. 



t ' 



t > 



and 



(d) 



B(l-9,9)J T _ 



(3) Finally, assume that ip'(0 + ) = and Ai^p(s) ^ for any s € R. Then, for any 5 < 1, 
we have the following intertwining relationship 



P-D 



j-jl, S pTifi^ 
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t > 



and 

Jt s , iP — B(1 — 6,S)J^. 

Proof. First, from the self-similarity property, we observe easily that the condition (|5.6j) is 
equivalent to the requirement that the kernel M^ B associated to the positive random variable 
J^g is injective. Since M.^ g {s) ^ for any s € R, we deduce that the multiplicative kernel 
M^ g is indeed injective, see e.g. [8j Theorem 4.8.4]. Next, note that ijj' e (0 + ) > and under the 
condition 5 < 9, 75 i eV ;/ (0 + ) > 0- Hence, from (|5.2j) . we deduce that, for any n > 1, 

The identity (|5.7p follows. Both processes being pssMps, the first intertwining relation follows 
from the criteria given above. The proof of the Theorem is completed by following similar lines 
of reasoning for the other claims. □ 

A nice consequence of the previous result is some interesting relationships between the eigen- 
functions of the semigroups of pssMps. Indeed, it is easily seen from the intertwining relation 
(15. ip that if a function / is an eigenfunction with eigenvalue 1 of the semigroup Pf then Af is an 
eigenfunction with eigenvalue 1 of the semigroup Qt- We proceed by introducing some notation 
taken from Patie |28j . Set ao(ifj) = 1 and define for non- negative integers n 

it 

k=l 

Next, we introduce the entire function which admits the series representation 

<x> 
n=0 

In [28\ Theorem 1], it is shown that 

(5.8) =l^(x), x > 0, 

where, for a smooth function /, the linear operator is the infinitesimal generator associated 

t 



to the semigroup Pf and takes the form 



Vf(x) = -xf"(x) + bf'(x)+x- 1 f(e' r x)- f{x) + xf'{x)rl { \ r \ <1} v{dr)-Kxf{x) 
1 Jo 



From the Feller property of the semigroup of X, we deduce readily that the identity (15 .8(1 is 
equivalent to 

e-*P t %(x) =X v ,(x), t,x>0, 
that is Xw, is 1-eigenfunction for Pf . Hence, we deduce from Theorem 15.21 the following. 
Corollary 5.3. (1) Let ip'{Q + ) < 0. Then, for any 5 < 9 + 1, we have the following identity 

D 9 '% 9 (x) =l T6 ^(x). 

(2) 7/V'(0 + ) < and 9 < 1, then 

(3) Finally, ifip'(0 + ) = 0, then, for any 5 < 1, we have 

D 1 '%(x) = l Ts ^(x). 
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We illustrate this last result by detailing some interesting instances of such relationships 
between some known special functions. 

Example 5.4 (Mittag-Leffler type functions). Let us consider, for any 1 < a < 2, the Laplace 
exponent ip(u) = (a(u — l/a)) a . We easily check that 6 = 1/a and we have tpi/ a {u) = (au) a . 
Observing that 

n im*o = i\^k) a = v{a( ? +i) \ n > i, 
k=i k=i 

and using the fact that the random variable J^ e is moment-determinate, we readily check, from 
the expression (|5.2p . that 

V(a(is + 1)) 



T(a)T(is + 1) 



The pole of the gamma function being the negative integers, the condition A4^p g (s) ^ is 
satisfied for any s£K. Moreover, we obtain 

where we recall that the Mittag-Leffler function £ a ,a is defined in the example 13.21 Next, for 
any 6 < 1 + 1/a, we have 

r(q)r(l/q + 1 - 5) ^ r(n + l/a + l) 

r M/aV>W r(l/Q + l) ^ r(n + 1/a + 1 - 5)Y{an + a) 

n=0 



2-^2 



^ (1, 1/a + 1), (1,1) 
(1, 1/a + 1 — 5), (a, a) 



where 2-^2 is the Wright hypergeometric function, see e.g. Braaksma [9j Chap. 12]. Hence, we 
have 

r( n \n e > s r M- v( ( 1 - 1 /« + 1 ).(M) \ 

Example 5.5. Now, for any 1 < a < 2, we set i)j{u) = u a and we note that ip'(0 + ) = 0. 
Proceeding as in the previous example, we get 

M^(s) = T a -\is + I) 

and hence the condition M^(s) 7^ is satisfied for any s € M. We have 



00 

v\ ) T a (n + 1 

fr,— n \ ' 



n=0 

and, for any 5 < 1, we write 



z r ^(x) = r(i - £ r rfr-M) . 



:o r(n + i-5)r«(n + i) 



Consequently, 

D 1 '%(x)=X T ,,oV( 2; )- 
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